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Kinetic equation for tachyons 
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The tachyonic regime of the quantum fluctuations of a self-interacting scalar field around its 
vacuum mean value is studied within a kinetic approach. We derive a quantum kinetic equation 
which determines the time evolution of the momentum distribution function of produced tachyonic 
modes and includes memory effects. The back-reaction of the quantum fluctuations on the vacuum 
mean field is taken into account, while their interaction is neglected. We show that the tachyonic 
modes do not correspond to real particles and contribute to the decay rate of the metastable vacuum 
state . 
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I. INTRODUCTION 


II. TACHYONIC REGIME 


Quantum fluctuations of a scalar field can enter a 
tachyonic regime where their frequency becomes imag¬ 
inary. Such regime can occur if the scalar field is either 
coupled to a strong stationary external potential |^,||,||j4) 
or strongly self-interacting with a potential exhibiting 
spontaneous symmetry breaking 110 - In the tachy¬ 
onic regime the system is essentially restructured. Its 
effective action develops an imaginary part |^, the fluc¬ 
tuations Hamiltonian becomes unbounded from below, 
while the Hilbert space of states acquires an indefinite 
metric |^. All these changes are indicative of a new, 
metastable phase. 

In the present work, we aim to derive a quantum ki¬ 
netic equation describing the production of the tachyonic 
modes for a self-interacting neutral massive scalar held. 
Particle production in the tachyonic regime has been ex¬ 
tensively studied so far in various models of spontaneous 
symmetry breaking |^. In these studies, the occupation 
number of produced particles has been estimated at the 
end of the metastable phase . Herein we suggest to study 
the full time evolution of the momentum distribution of 
the tachyonic modes using a kinetic description. 

The decay of the metastable vacuum state has been 
discussed in different ways, including the semiclassical 
approach |^, the classical lattice field theory |1I|, the 
two-particle irreducible effective action formalism [^ . 
Our approach is based on the canonical quantization of 
the tachyonic modes. 

The plan of the paper is as follows. In SecJ^ we in¬ 
troduce the model and identify the tachyonic regime. 
In Sec. we perform the quantization of the tachyonic 
modes. A quantum kinetic equation is derived in Secp^. 
We conclude with summary in Sec.^ 
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We consider a general scalar field model with the La- 
grangian density 

^ ( 1 ) 

where V{(p) is a self-interaction potential which contains 
orders and higher without derivative terms and m is 
the mass of the scalar field. The model is defined in a 
finite volume L^, —L/2 < Xi < T/2, i = 1,2,3. The 
continuum limit is f 

From ( 0 ) we obtain the Klein-Gordon type equation of 
motion for the held p(x,t): 

{a+ Tn'^)p = J = (2) 
op 

where the non-linear current J is also determined by the 
self-interaction. 

Following the mean-field approximation, we decom¬ 
pose p{x, t) into its space-homogeneous vacuum mean 
value (fit) = {p{x, t)) and fluctuations x 

p{x,t) = (j)it) + xix,t) (3) 

with ixix, t)) = 0. The mean field is treated as a classical 
background field defined with respect to the in-vacuum 
|0) as 


(fit) = ipix,t)) = ^ J d^xi0\pix,t)\0), (4) 

so that in the limit t —f —cx) (fit) —*■ 0, while the fluctu¬ 
ations are quantized and take place at all times. 

Using Eq.(^ provides the following decomposition for 
the current 

Ji^ + X) = Ji<t>) + ^^^^X + ^i<t^,x), (5) 

where x) includes terms of second and higher orders 

in X, 
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>/(</>, x) 


152j(0) 2 

2^^^ + 


( 6 ) 


Substituting Eq.(^) also into Eq.(||) and taking the mean 
value (...) yields the vacuum mean field equation 


ij) + (j) - J{(j)) = {J), (7) 


where the overdot indicates the derivative with respect 
to time, while the equation of motion for the quantum 
fluctuations reads 


the time dependence of the vacuum mean field 4>. If (j) 
oscillates, then the same momentum mode can change its 
nature during the time evolution. 

Eqs. (I) and & are exact, self-consistently coupled 
and include back-reactions. The vacuum mean held mod- 
ihes the equation for huctuations via a time dependent 
frequency, while the huctuations react back on the vac¬ 
uum mean held via the source term (J) in Eq. ( 0 ) and on 
the huctuations themselves via the “external force” term 
F^{k,t) in Eq.(|^. 


(□ -I- — J — (J) 


with 


‘e// 


5(tP ' 


( 8 ) 

( 9 ) 


r2 T/" 

Eor > 0, the effective mass squared is positive at all 
times. However, if < 0, becomes negative for 
> m? indicating a tachyonic regime. 

In terms of the Eourier components x(fc, f), Eq.(||) 
takes the form 


X(fc, t) + Lol{t)x(k, t) = F^ik, t), (10) 

where 

Fxik,t) = J{k,t) - VV{J)S^^g (11) 

and J{k, t) is the Eourier transform of the current J{x, t), 

J{k,t) = j (12) 

while 

ujlit) ^k^^ (13) 

is the time-dependent frequency squared of the huctua¬ 
tions . In the tachyonic regime, w^(t) can be negative. 

Whether the system evolves in the tachyonic or non- 
tachyonic regime is dynamically hxed by the time- 
dependent critical momentum: 


III. QUANTIZATION 

The Hamiltonian density corresponding to (||) is 

77 = -h i(V(/5)2-f + ^(V^), (15) 

where tt is the momentum canonically conjugate to tp. 
With the decomposition for the potential 

Vi(l) + x) = Vicj}) - J{4')x + ]^{mlff-m^)x^ 

+ n</>,x), (16) 

orders x^ and higher being included into V{(j), x) i we de¬ 
duce from (H) the Hamiltonian density governing the 
dynamics of the huctuations 

+ V((/>,x)- (17) 

In terms of the Eourier components x(k, t) and 7rx(k, t), 

the huctuations Hamiltonian reads 

Hx= j dPxHx 

= ^ XI (4(^’^)’^x(^>^)+^^fc(0x^(fc,i)x(^,i)) 

+ i {^lik,t)'^xik,t) - i^kit)xHk,t)xik,t)j 

k 

+L^/W{k = 0,t), (18) 




0 , 


— TO 


otherwise 


(14) 


All momentum modes below k^ are tachyonic. Eor 
> —m?, the critical momentum is zero, since the 
frequency is always positive and no tachyonic modes can 
appear. 

The system can enter the tachyonic regime in diherent 
ways: gradually when the critical momentum becomes 
nonzero very slowly, or discontinuously when tachyonic 
modes appear suddenly [|^,|^,^ on a short time scale. 
In any case, the critical momentum changes in tune with 


where = —io\ > 0 for fc ^ and 


x\k,t) =x{-k,t), (19) 

T^x{k,t) = TTx{-k,t) ( 20 ) 

for all momentum modes. 

The non-tachyonic and tachyonic modes contribution 
to the Hamiltonian (0) represents a collection of positive 
and inverted (repulsive) oscillators, respectively. Both 
types of modes are coupled. Their interaction is de¬ 
scribed by the last term in Eq.([l^), V{k,t) being the 
Eourier transform of the potential V{x,t). 
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( 29 ) 


For the standard, non-tachyonic modes, we introduce 
the annihilation and creation operators by 

Xik,t) = rj(t)a(fc,t) +Ti{t)a\-k,t), (21) 

and 


TT^{k,t) = -iujk(t) Tj:(t)a{-k,t)-T*M)a\k,t) , (22) 


where 




1 

\/2wfe(t) 


exp{-z0fe(wfc,t)}, 


(23) 


and 0fc(wfe,t) is a phase which in the in-limit takes the 
form uj it = \/k'^ + m?t. 

Eqs.(pT|) and (^^ are well-known expressions for the 
real frequency oscillations. The first term in the Hamil¬ 
tonian ([l^ ) - we denote it by - becomes up to a 
c-number: 

^x‘= E (24) 

k 2>fc^2 


where N'^*{k,t) = a^k,t)a{k,t) is the non-tachyonic 
modes number density operator. 

For the modes with k ^ < k^, LOj^ = ±Wfc = 


— k"^ and one of the phase factors in the ansatz 
(^), r^(t) or rt(t), grows exponentially in time. In¬ 
stead of oscillations we have an exponential growth of 
long wavelength quantum fluctuations with momenta 
k ^ < k^. This is the so-called tachyonic instability 


Making the transition ujk —> Vk in Eq.(22) yields the 
following ansatz for the negative frequency squared fluc¬ 
tuations 


X{k,t) Xt{k,t) 


- ^^(e’^'=a(fc,t)-ke fc,t)), 

(25) 

where ‘dk{vk,'t) = —*0fc(wfc,t). Introducing 


tTi(fc, t) = / cosh-dfe • {a{k, t) + a\—k, t)] 

y/2vk V 7 

, (26) 

a 2 {k,t) = — — ^i=sinh'dfe • (^a{k,t) — a\—k,t) 

), (27) 

which obey the hermiticity condition <7^2) 
CTi( 2 )(— A, t), we rewrite Eq.(p5[) as 

(fc, t) = 

Xt{k,t) = cri(fc,t) -1- i(T2{k,t) 

(28) 


with xj(fc,t) ^ xti—k,t), i.e. the ansatz ( p^ ) is non- 
Hermitian. 

The canonically conjugate momentum is transformed 
as 


Trx{k,t) TT^^t{k,t) = TTcn{k,t) + iT:^2{k,t), 
where 


7ro-i(fc, f) = zzfc cothdfc • (30) 

7ro.2(fc,f) =-r'/ctanh-dfe • cr|(fc,t). (31) 

The commutation relations for ai ,cr2-fields are 


ai{k,t),a2{p,t) 


2 

-— sinh 2dfc • 
2l2k 




(32) 


all other commutators vanishing. 

Analytically continuing the ansatz (ID in the fre¬ 
quency to imaginary values leads therefore to a non- 
Hermitian field. This is not acceptable because we re¬ 
quire the hermiticity conditions, Eqs.(19)-(20), to be 
valid for all momentum modes and at all steps of our 
consideration. In addition, such non-Hermitian field is 
known to violate causality . 

To define the Hermitian tachyonic fluctuations we can 
use either cri(fc, t) or a 2 {k, t) instead of Xt{k, t). Without 
loss of generality, we choose ai{k,t) and introduce the 
field 


at{k,t) = 


1 


-ai{k,t). 


cosh dfc 

Its canonically conjugate momentum is 


(33) 


Tr^^t{k,t) = -7r<^,i(fc,f). (34) 

cosh x/k 

With Eqs.(p^) and (|3^), the second term in the Hamil¬ 
tonian (p8|) takes the form 


Hx= Y. (35) 

k 


where 

N*{k, t) = — — (^a^{k, t)a\—k, t) + a{—k, t)a{k, t)^ . (36) 


Since the spectrum of an inverted oscillator is purely 
continuous, the tachyonic modes are not really “parti¬ 
cle” ones [Q. In contrast with the case of the standard, 
non-tachyonic modes where the eigenfunctions of co¬ 
incide with those of the number operator, the tachyonic 
modes are not eigenoperators of 

JV*= Y (37) 

k 


namely 


N*, a{k, t) = a\—k, t), 
N*^,a\k,t) =—a{—k,t), 


(38) 

(39) 
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so that a{k,t),a'^{k,t) in Eq.(^ can not be viewed as 
creation and annihilation operators. 

However, once complex values are allowed for energy, 
the particle interpretation can be kept for the tachyonic 
modes as well . Let us introduce 

ai(2)ik,t) = + ^-^a{k,t), (40) 

where the upper signs correspond to the subscript 1 and 
the lower ones to 2. These new mode operators are Her- 
mitian and fulfill the algebra 


ai{k,t),al{p,t) 


ai{k,t),al{p,t) 


a2(fc, t), a\{k, t) 


= 0 , 


= i5: 


k,p' 


(41) 

(42) 


The Fock representation for the algebra (41)-(42) is con¬ 
structed by using an indefinite metric . Indeed, if |0; t) is 
an instantaneous vacuum state defined by 


N^\na2; ki, fc2, k„; t) = n\na2] ki,k2, t), (51) 

i.e. plays the role of the “number operator”. 

In the space with indefinite metric, the Hamiltonian 
iL* is pseudoadjoint ^ and its eigenvalues are imaginary. 
If |e; t) is an eigenstate of with eigenvalue e, then for 
the state al{k,t)\e;t) we obtain 

H^al{k,t)\e;t) = (^e + iL'k^al{k,t)\e;t), (52) 

i.e. al{k,t)\e;t) is also an eigenstate of iJ* with the 
eigenvalue shifted by Wk- Neglecting for a moment the 
third term in the right-hand sise of Eq.(p^, we see that 
the eigenvalues of the total Hamiltonian are 

complex, the corresponding eigenfunctions representing 
unstable states. 


IV. KINETIC EQUATION 


ai{k,t)\0;t) = a2{k,t)\0-,t) = 0 for (fc^) > 0, (43) 
where {ki) = (fci, A: 2 , fca), then for the excited states 

|ai(2);fc,t) = al(2)(fc,i)|0;t), (fci) > 0, (44) 

the inner product is vanishing or imaginary, 

(ai(2);fc,t|ai(2);p,t) = 0, (45) 

{ai;k,t\a 2 -.,p,t) = (46) 

The indefinite inner product is related to the existence 
of associated eigenvectors of |Q. 

The density of N* becomes 

N\k,t) = -iNl{k,t) 

= -^(«i(^>^)«2(fc,i) + al{k,t)ai{k,t)'^, (47) 


In the mean-field approximation, the quantum fluctua¬ 
tions are treated perturbatively. This is necessary, in par¬ 
ticular, for the derivation of the kinetic equation. One of 
basic points of the kinetic formulation is the postulate of 
asymptotic completeness 1^ . The postulate specifies the 
set of possible states of the system in the infinite past as 
a complete set of states of freely moving non-interacting 
particles. For the system with a self-interaction, this pos¬ 
tulate can be applied only in a few lower orders of per¬ 
turbations when the interaction potential vanishes in the 
in-limit due to the vanishing of the vacuum mean field . 
In higher orders, the quantum fluctuations dominate, the 
corresponding terms in the interaction potential surviv¬ 
ing in the infinite past. 

We limit our consideration to the third order term in 
V{(j), x) neglecting all higher orders. In addition, we use a 
Hartree-type approximation that in the second and third 
orders consists of the factorization 


cti{ 2 ){k,t) being eigenoperators of 


Y. 2iV‘(fc,t) 

k 

( fci )>0 


with real eigenvalues, 

Nl,a^2){k,t) 


Toti(2){k,t). 


(48) 


(49) 


For the instantaneous vacuum, N^\0; t) = 0, while for the 
excited states counts excitations. For the state 


1710 ^ 2 , k \, k2 ,..., kji , t) 

= a\{ki,t)a\{k2,t) ■ ... ■ a|(fc„, t)|0;<), 

all (^1,2, k^i^i) 0 , 


^ (x^), X^ ^ 3(x^)x. (53) 

For the non-tachyonic modes, the form of the kinetic 
equation is well-known and was given in different models 
It determines the time evolution of the occupa¬ 
tion number density 

AA”*(fc,t) = (0|A"*(fc,t)|0) (54) 

which defines the number of particles of a given state 
characterized by the momentum fc ^ > k^^ at time t. An 
increase in the occupation number density is interpreted 
as particle production. 

Herein we focus on the time evolution of 

M\k,t) = {0\Nl{lt)\0) (55) 


for instance. 


( 50 ) 
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which defines the momentum distribution of the tachy¬ 
onic modes. We start with the tachyonic Hamiltonian 
equations of motion 












at{k,t) = Trlt{k,t), (56) 

TT^^tik, t) = vial (k, t) + J{-k, t), (57) 


where the current J(—k,t) represents the self-interaction 
potential contribution. With the factorization the 

self-interaction potential and current take the form 

V(k = 0,t) =-^^-^^{X^)at(0,t) (58) 

and 

= (59) 

respectively. Using the relations 

Oii{2)(k,t) = ]J^{^t(k,t) T ^Trl^tik,t)'^ (60) 

yields then the equations for ai( 2 )(fc,t): 


ai(2)(fc,t) ± Vkai(^2){k,t) - -^a2(i){k,t) 


= T 


J(k,t). 




(61) 


Taking next the time derivative of J\r*(k,t) we find: 


N\k,t) = ^^{c^(k,t) + C2(k,t)), (62) 

where we have defined the time-dependent one-particle 
correlation functions 

C'i(2)(fc,f) = (0|ai(2)(-fc,t)ai(2)(^,i)|0)- (63) 

Since (x(a:, t)) = 0, the vacuum expectation values for the 
zero momentum mode operators (Tt(0, t) and 7rCT_t(0, t) are 
equal to zero, and , as a result, the current J(k,t) drops 
out of Eq. ( |6^ . 

The functions Ci( 2 )(k,t) obey the equations 

Ci{2)(k,t) = —A/'*(fc,t) T 2z/feCi(2)(fc,t). (64) 

Their formal solution is 

Cl(2)(fc,t)= f dt''^N\k,t') 

Jto Mt') 

X exp{±2«(t') - df(t))}, (65) 

where 

= [ dt’vk(t') ( 66 ) 

Jto 

and to is a moment of time at which the tachyonic regime 
starts. If to = —oo , then the in-vacuum can be chosen 
as an initial state of the system. 


Although we have not assumed that the frequency Vk 
varies adiabatically slowly in time and the phase dk in 
the ansatz ( ^) is a general function of cok and t, it is just 
the “adiabatic” phase (|^, i.e. the phase which looks 
exactly like the one in the adiabatic case, that enters 
this solution. Substituting it into Eq.(|^), we obtain a 
closed equation for 


— r 

Vk Jto 

X cosh hdfit') - 2df(t) 




(67) 


This kinetic equation determines the time evolution of 
the momentum distribution of the tachyonic modes pro¬ 
duced in the fluctuations of the scalar field. As seen from 
the definition (^), the tachyonic modes production is 
symmetric in the momentum space, Af*(k, t) = Af*(—k, t) 
for all times t. 

Eq.(^^ has non-Markovian character due to the ex¬ 
plicit dependence of its right-hand side - the source 
term — on the time evolution of J\f*'(k,t) and therefore 
involves memory effects starting from to- For the real 
particle modes, the source term is known to contain a 
time integration over the statistical factor (1 ± 2Af(k, t)), 
where the plus sign corresponds to bosons and the minus 
one to fermions |^. For the tachyonic modes , this fac¬ 
tor reduces to 2J\f(k, t) reflecting once more the fact that 
tachyons are not real particles . 

In our approximation, the vacuum mean field equation 
becomes 

(f, + rn^(l) = + (68) 

where the (x^)-term represents the back-reaction of the 
fluctuations on the vacuum mean field and provides 
damping of the oscillations of (p. The initial conditions 
for both Eqs.d^) and (^) are specified by the model 
under study. 

The vacuum mean value of is given by 
(X^) = (^\xik,t)x{-k,t)\0) 

+ ^ X! {0\crt(k,t)at(-k,t)\0), (69) 

k 

the bi-linear operator expressions here being assumed to 
be normal-ordered with respect to the instantaneous vac¬ 
uum state. Both types of modes, tachyonic and non- 
tachyonic, contribute to (x^), so a proper inclusion of 
the back-reactions effects can be achieved only by the 
complete treatment of all momentum modes. 

Taking the vacuum expectation value of the fluctua¬ 
tions Hamiltonian yields 
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(70) 









where the non-tachyonic modes contribute to the energy 
of the metastable vacuum state 

E^= Y. (71) 

while the tachyonic ones to its decay rate, 

T^ = 2 Y (72) 

k 


V. SUMMARY 

For the model of a self-interacting scalar field, we have 
derived a non-Markovian quantum kinetic equation de¬ 
termining the momentum distribution of the tachyonic 
modes. These modes are produced in quantum fluctu¬ 
ations of the scalar field around its vacuum mean value 
when the system is in a metastable phase. The kinetic 
and vacuum mean field equations are coupled, the lat¬ 
ter including the back-reaction term, while the collisions 
effects are neglected. 

Despite the fact that the fluctuations Hamiltonian is 
not bounded from below in the tachyonic regime, the con¬ 
servation of energy prevents any catastrophic production 
of tachyons. If the system starts in a false, metastable 
vacuum state and then undergoes the transition to a 
lower energy density , stable one, the tachyonic regime 
stops as soon as all the potential energy of the false vac¬ 
uum state is transferred into the quantum fluctuations. 
We have shown that the tachyonic modes contribute to 
the decay rate of this state, so their intensive production 
results in its rapid decay. 

The kinetic equation obtained is hoped to be useful 
for the numerical study of the tachyonic modes produc¬ 
tion in various problems, in particular, in cosmology ||] 
and heavy-ion collisions |p^ , p^Jl^ . The complete study 
requires the inclusion of higher orders effects when the 
quantum fluctuations interact with each other. Its real¬ 
ization within the kinetic formulation would provide fur¬ 
ther insight into the dynamics of the tachyonic regime. 
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